Peierls distorted chain as a quantum data bus for quantum state transfer 
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We systematically study the transfer of quantum state of electron spin as the flying qubit along a 
half-flUed Peierls distorted tight-binding chain described by the Su-Schrieffer-Heeger (SSH) model, 
which behaves as a quantum data bus. This enables a novel physical mechanism for quantum 
communication with always-on interaction: the effective hopping of the spin carrier between sites A 
and B connected to two sites in this SSH chain can be induced by the quasi-excitations of the SSH 
model. As we prove, it is the Peierls energy gap of the SSH quasi-excitations that plays a crucial 
role to protect the robustness of the quantum state transfer process. Moreover, our observation also 
indicates that such a scheme can also be employed to explore the intrinsic property of the quantum 
system. 
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I. INTRODUCTION 



In many protocols of quantum information processing, 
it is crucial to transmit the quantum state of qubits with 
high fidelity T] . While various schemes of quantum state 
transfer (QST) were proposed and demonstrated exper- 
imentally J3] for the optical system even with atom en- 
semble people have tried the best to implement this 
task based on the solid-state systems, which are believed 
as the best candidates for the scalable quantum comput- 
"2 01 • It was recognized that the quantum spin chain 
and the Bloch electron system j^] with the artificial 
nearest neighbor (NN) couplings can be used as quantum 
data buses to transfer quantum information perfectly. In 
Ref. Q, we discovered that the spectrum-parity match- 
ing is responsible for the perfectness of most protocols of 
QST. 

However, these solid-state based schemes for imple- 
menting quantum data bus are too artificial with very 
specially designed NN couplings and only single-particle 
cases were considered. Though it is believed that such an 
engineered quantum system will be realized in the future 
experiments, at least, this difficulty should be overcome 
in theoretical aspect. To this end, we have tried to em- 
ploy the higher-dimensional system or the complicated 
quantum network, such as the spin ladder jlOl |. Though 
the spin ladder is not an ideal medium serving as a perfect 
data bus, it shows that the long-distance QST is possible 
via the multi-exciton system due to the existence of the 
spin gap, which is much larger than the coupling strength 
connecting two separated qubits to the two ends of the 
ladder respectively. 
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FIG. 1: (a) The schematic illustration of the system. Two 
sites A and B connect to the two sites of PDTBC with the 
distance I being (b) odd and (c) even. 



A medium can be a robust quantum data bus if the fi- 
delity of transferring a quantum state between two sites 
attached to it approaches unity at a finite temperature. 
These consideration motivates us to use more natural 
energy gapped materials to act as a quantum data bus. 
In this paper, we study the function of quantum data 
bus for the conducting polymers (polyacetylene), mod- 
eled as a Peierls distorted tight-binding chain (PD TBC) 
or the Su-Schrieffer-Heeger (SSH) model [H It is 
well known that, taking the coupling between electron 
and photon into account, such a system exhibits an in- 
stability against the lattice distortion, which induces an 
energy gap for the spectrum in half-filled case as a re- 
sult of dimerization. Therefore, the SSH chain is a good 
candidate of quantum data bus for fermion, which is just 
similar to the spin ladder for spin 10]. Furthermore, such 
a practical system allows us to investigate the case with 
larger particle number by performing analytical analysis 
for the two distant separated probe points AB since the 
SSH model can be solved exactly. 



2 



As a protocol for the quantum communication with 
always-on interaction, two points A and B are connected 
to two sites of the SSH chain (see the Fig. 1(a)). For 
the half-filled case, when the connections of A and B 
with SSH chain are switched off, the ground states are 
two-fold degenerate due to the vanishing on-site chemical 
potential. When the connections switch on, the two de- 
generate levels split while the eigenstates remain nearly 
half-filled for the subsystem AB due to the energy gap 
of the SSH chain. It leads to the effective interaction 
between A and B. Thus the level spacing of two lowest 
eigenstates directly corresponds to the effective hopping 
integral between A and B and determines the validity of 
quantum state transfer via such a medium. Generally, if 
the energy gap is large enough, an effective Hamiltonian 
Hab of these two sites is induced by this quantum data 
bus to perform perfect QST at low temperature. Thus 
the validity of Hab and the behavior of the effective hop- 
ping integral (level spacing of two lowest eigenstates) as 
the distance between A and B increases are crucial for the 
quality of a quantum data bus. To demonstrate the prop- 
erties of our scheme, we compare the eigenstates of Hab 
with those reduced density matrices from the ground and 
the first excited states of the total system. We find that 
the Peierls energy gap of the SSH chain plays an impor- 
tant role in protecting the robustness of the QST. The 
main conclusions are achieved both with the analytical 
study and numerical calculation. 

This paper is organized as follows. In Sec. II, the 
model setup of our protocol and the spectrum of SSH 
chain are introduced. As a reasonable approximation up 
to second order, the effective Hamiltonian Hab with re- 
spect to the two separated points A and B is deduced by 
using the Frohlich transformation [13 to "remove" 
the degree of freedom of SSH model in Sec. III. In Sec. 
IV, the density matrices of qubits A and B for the ground 
and the first excited states of the whole system are calcu- 
lated to demonstrate the quantum entanglement of them 
over a long distance. The QST scheme is also investi- 
gated numerically in the region of the crossover between 
two types of dimerization. It shows that, at this transi- 
tion point, the QST becomes fastest, but with the similar 
lost of fidelity to what has been discussed in the Ref. [23| . 
In Sec.V, we discuss the quantum decoherence problem 
due to the quasi-excitations of the SSH chain at finite 
temperature. The summary and remarks are presented 
finally in Sec. VI. 



II. MODEL SETUP OF QUANTUM DATA BUS 
BASED ON THE SSH CHAIN 



In a polyacetylene, cr-electrons, which are localized be- 
tween the two bonded nuclei, connect one C atom with 
the neighboring two C and one H atoms to form an ele- 
mentary linear configuration; 7r-electrons, which are less 
localized, behave as Bloch electrons. Such a system can 
be treated by using tight-binding approximation and tt- 



electrons can be the candidate as a carrier of quantum 
information ||9j. Here, we also consider the coupling be- 
tween electrons and dispersionless phonons represented 
by local oscillators. 

Let Un be the displacement of the nth CH (the unit 
of a polyacetylene) from the equilibrium position. Then, 
the tight-binding Hamiltonian reads as 



H 



SSH 



h.c. 



He 



where 



H 



c 



E 



(1) 



(2) 



describes the lattice motion. Here, = .90 ^ 

A (un+i — Un) is the hopping integral and A is the 
electron-lattice coupling constant. n„ denotes the canon- 
ical momentum operator conjugate to u„ and M is the 
effective mass of the unit of a polyacetylene (CH). K is 
the spring constant due to cr-electrons, and Cj^a (c] '^^ 
annihilation (creation) operator of 7r-electrons on site j 
with spin cr = ±1. This Hamiltonian is proposed by Su, 
Schricffer and Heeger to describe the dimerization 
phenomenon of the polyacetylene in association with its 
conductivity. 

In our scheme, we consider the QST between sites A 
and B connected to the two sites of the SSH chain. Our 
purpose is to transfer a qubit state, which is in a super- 
position state of electron spin up and spin down, through 
the hopping of an polarized electron from the site Aio B 
(see the Fig. 1(b)). Since the Hamiltonian does not 
contain spin dependent interaction, the spin polarization 
of every electron is conservative. In this sense, the spatial 
motion of electrons along the chain will carry a spin state 
from one location to another. Thus, we can study the fi- 
delity of spin state transfer by considering the fidelity of 
charge state transfer. The connection of sites A and B 
to the SSH chain is described by the Hamiltonian 



Hj 



(3) 



where CA,a and CB,a (c^ „ and ^) are annihilation (cre- 
ation) operators of electrons on sites A and B with spin 
cr, while Iq and Iq + I denote the connecting sites in the 
SSH chain. 

In Fig. 2, we intuitively demonstrate how the above 
Hamiltonian can lead to the effective hopping of elec- 
tron between the additional sites A and B. We will also 
show how the effective hopping integral depends on the 
distance between them. The original energy degeneracy 
at = of the subsystem formed by A and B will be 
removed by the switching on g and then split into two 
sub- levels with the level spacing 2geff. Here, geff is the 
effective hopping integral, which can be obtained analyt- 
ically as follows. 
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FIG. 2: (a) The single particle energy spectrum of dimerized 
polyacetylene. (b) The eigenstates of the effective Hamilto- 
nian about points A and B. 



To consider how a SSH chain can be used as a quan- 
tum data bus to transfer a quantum state, we first sum- 
marize the known results about diagonalization of the 
SSH Hamiltonian. The most important discovery is the 
prediction about the Peierls transition happens when 
the polyacetylene is dimerized {-!)"-' uo. It 

minimizes the energy of the one-dimensional electronic 
gas when the lattice vibrates slowly so that the Born- 
Oppenheirmer approximation can determine an instan- 
taneous eigen energy of electron for a fixed lattice config- 
uration. After dimerization, the Hamiltonian then reads 



HssH - - ^ .90 [1 - (-1)" S] (c,\+i,,c„,, + /i.e.) + Eo, 

71,(7 

(4) 

where S = 2Xuo / go denotes the distortion of the hopping 
integral and Eq is a constant. Taking periodic boundary 
condition, it can be diagonalized as 

HsSH = X! ^'^ {^k,a^k,a " ,„/3fc,o-) , (5) 



with the dispersion relation 



efc = 2go\/ cos2 ^ + sin^ ^ 



(6) 



for the excitations described by the fermion operators 

N/2 



and 



N/2 



where N is the number of CH in polyacetylene chain 
(or the length of the chain), k — Attiti/N, m = 
0,1,2, ...,7V/2-l, and 

e''^" =^^[{l + 5) + {l-5)e-'^]. (9) 

The single-particle spectrum © is illustrated in the Fig. 
2(a) and the energy gap between the two bands is 



2A = 2min{efc} = Ag^S. 



(10) 



In next section, we will show that such an energy gap is 
necessarily required for a desirable robust quantum data 
bus for QST. Usually, as illustrated in Fig. 2(b), the 
large energy gap above the lowest two levels can provide 
a kind of "quantum protect" for the QST via the vir- 
tual excitations of the quantum data bus, which induce 
a long range hopping of electron between sites A and B in 
the low temperature. The dense continuity of spectrum 
above the gap ensures that the strength of the effective 
hopping integral across A and B is so strong that the fast 
entanglement can be generated. In the following, we will 
demonstrate that such a dimerized gap play s a same role 
for QST as the spin gap in spin ladder 



play s 



III. EFFECTIVE LONG RANGE HOPPING 
INDUCED BY VIRTUAL QUASI-EXCITATIONS 

To deduce the effective Hamiltonian about the indirect 
coupling between two attached sites A and -B, we utilize 
the Frohlich transformation, whose original approach was 
used successfully for the BCS theory of superconductiv- 
ity. The effective Hamiltonian -ffe// = UHU~^ can be 
achieved by a unitary operator U — exp {—S), where H 
is the Hamiltonian of the total system with the pertur- 
bation decomposition 



SSH 



(11) 



In the second order perturbation theory, we require the 
anti-Hermitian operator S obeys 

Hi + [HssH,S] = 0. (12) 

Thus the effective Hamiltonian can be approximated as 

1 



^e// = HssH + 2 [Hi, S] 



(13) 



There are two ways to connect the two sites A and B 
to the SSH chain, with (without) the mirror inversion 
symmetry as showed in Fig. 1(b) (Fig. 1(c)). Here, 
only the case with both Iq and I being odd is discussed. 
The result can be used to all cases, because the system 
is invariant under exchanging A and B. As the solution 
of Eq. H12|l , S can be expressed explicitly as 

S - -^Ef (14) 



N "-^ ek 

k.a 

,fc!o+l_,0 t 



[ak.cr + (3k.a) - h.c. , 



which leads to the effective Hamiltonian 
Heff = Hab+Ho 



with 



Here 



[Hab,Ho]=0. 



Has = XI 9eff {c\ aCB,a + C^B.a'^A,a 



(15) 
(16) 
(17) 



denotes the effective hopping of electron between A and 
B with the strength 



9eff 



k 



(18) 



2N 



while 

k.(7 

X X — U (fc, fc') (4,o-"fc','^' - Pl,aPk',a') 
k,a,k' ,cr' 

+ B (fc, fc') iPl^ak'.a' - aUf3k'.a') + /i.e.] (19) 
describes the dynamics of the data bus, where 



'o + i 



1 _ e*(«fc-efc') 



(20) 



Straightforward calculation shows that, in the thermo- 
dynamic limit, N — > oo, the hopping constant becomes 



9eff 



go 1 + S \1 + S 



(21) 



Here, we have used the one-dimensional integral to re- 
place the discrete sum as (2/-/V) — > J^^ dk/2'K. 

It needs to be pointed out that the effective Hamilto- 
nian Hab works well only in the half-filled case. When 
the lower quasi-band is filled fully, the next excitation 
should jump over the gap and then virtual excitations 
can effectively induce the effective coupling. To demon- 
strate the above result about the effective long range hop- 
ping induced by the SSH chain, we schematically sketch 
the virtual process in Fig. 3, where the quanta of data 
bus between the two sites are exchanged. If the energy 
gap of the PDTBC is much larger than the connection 
couplings, A 3> g, the virtual transition is essentially 
second order process. There are two kinds of virtual pro- 
cess for an initially half- filled state of the PDTBC. The 
first one happens via the exchange of the electron in the 
upper quasi-band, while the second one is based on the 
exchange of the hole in the lower quasi-band. We would 



M 



(a) 




(b) 



FIG. 3: Schematic illustration for the two processes of the 
second order perturbation. The virtual transitions of (a) elec- 
trons between A{B) and the upper quasi-band and (b) holes 
between lower quasi-band and A{B). These processes lead to 
effective hopping of electron between A and B directly. 



like to remind that the two quasi-bands will become con- 
tinuous in the thermodynamic limit. Nevertheless, there 
still exists indirect interaction between A and B as the 
effective results of the second order process. 

Now we consider the scheme using the SSH chain to 
transfer the quantum state via the above effective virtual 
excitation process. Let Alice hold an electron with spin 
state 



|^)=cos-U)+e'^sin-|T) 



(22) 



at the site A, where jt) (|J,)) denotes the spin up (down) 
state. Thus, the initial state (0)) = |<y5)^ <8) |0)^ of the 
total system is 



1* (0)) = (^cos -c\^^ + e'* sin -c\^^j 10)^^ . 

Here, |0)^ (|0)^^) denotes the empty state, i.e., there is 
no electron at site A (both sites A and B). At the instant 



2|ffe//|' 



(23) 



the total system evolves into a new factorized state 

I*M) = |o)^®|<p)b ot 

1* (r)) = (^cos ^-cl^^ + e'* sin ^-c^^^^ \0)^g , 

to realize a perfect quantum swapping. Then Bob at 
the site B can receive an electron with spin state \(p). 
Here, t determines the characteristic time of quantum 
state transfer between two locations A and B, so the 
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FIG. 4: The magnitude of effective hopping integral g^ff for 
the finite systems obtained by numerical exact diagonalization 
(solid dot) and the approximate analytical results from Eq. 
1211 (line). Here only odd I is plotted. It shows that geff 
decays not so fast in this range of I and the exact numerical 
and approximate analytical results are in agreement with each 
other very well. 



tonian Hab by comparing the eigenstates of Hab with 
the density matrix reduced from the ground and first ex- 
cited states of the total system Notice that the 
Hamiltonian does not contain the spin-spin interac- 
tion term. Then we can regard the system as a spinless 
fermion system and then apply the result obtained feasi- 
bly for the spinless case to the original system. Therefore, 
in the following discussion we will ignore the spin degree 
of freedom for simplicity. 

We define the states \n^,n+) by 



|0,1)ab 



|o) 

|0> 



AB ' 



AB 



(24) 



to denote the eigenstates of Hab in the half-filled sub- 
space, where a^, are spinless fermion operator at sites 
A and B. Here, n_ (n+) is the number of particles in 
the anti-bonding (bonding) state. Then we calculate the 
"fidelity" of eigenstate defined by P„_„^ = Tr{pfipn_n+) 
or 



(25) 



behavior of |ge//| vs. / is crucial for quantum information 
transfer. 

It is desirable that the above scheme for QST based on 
the SSH virtual excitation can work well over a longer dis- 
tance, namely, the effective hopping integral ffe// should 
not decay too fast as the distance of two points AB in- 
creases. According to Eq. if^ . \geff\ decays exponen- 
tially as I increases. However, ^{1 — 6) / {1 + 6) closes 
to 1 for small S. So \geff\ does not decay rapidly for any 
finite To demonstrate this, numerical simulation is 
performed. The value of |ge//| is computed from the en- 
ergy spacing between the ground and first excited states 
of the total system by exact diagonalization method for 
the finite system with N = 500, I ^ lOn + 9, n e [0, 7], 
5 = 0.01, g = 0.01, and go = 1.0. The exact numerical 
result and the analytical result obtained from Eq. H21(l 
are plotted in Fig. 4. It shows that they are in agreement 
with each other well and |5e//| exhibits power law decay 
for finite I. It indicates that the characteristic time r is 
proportional to the distance I, which is crucial for the 
scalable quantum information processing. 



IV. REDUCED DENSITY MATRIX AND 
ENTANGLEMENT 

In the above section, we studied the behavior of |5e//| 
changing as the distance / increases. It determines the 
speed of the QST via this SSH chain. However, all the 
conclusions obtained should be based on the fact that the 
effective Hamiltonian Hab given by H17|l is valid. In this 
section, we investigate the validity of the effective Hamil- 



Obviously, it can be described as the expecta- 
tion value of the effective density matrix p„_n+ = 

(n_, n+l in the exact eigenstates 
of the total system with the reduced density matrix 
Pfi = TrssH{\'n-,n+) {n^,n+\). Here, \n^,n+,T]) = 
|n_,rt+)^g l'?)ssH denotes the eigenstates of ((TT|l with 
g = 0, while \'r])ssH denotes the groundstate (77 = 0) and 
excited eigenstates {ij — 1,2,...) of the half- filled SSH 
chain. State |n_, is the ground or first excited states 
of the whole system (the Hamiltonian Hll|) with 5 7^ 0), 
which are spanned by the states possessing the same par- 
ity as the states |n_, n+; 0) — |n_, (g) |0)gg^. 

On the other hand, it is known that the eigenstates 
of Hab with unity concurrence are maximally entangled 
states Thus by calculating the mode concur- 

rences of two sites A and B for the ground and first 
excited states of the total system Ullf) . one can verify 
whether the exact solutions can be well approximated 
by that from the effective Hamiltonian. If "fidelities" or 
concurrences is near unity, it can be concluded that the 
effective Hamiltonian Hab is indeed valid for the perfect 
QST. 

In the following discussion, we try to explicitly express 
the exact eigenstates of the total system with 

the help of the Gellmann-Low theorem • In quantum 
field theory, this theorem is used successfully to obtain 
the ground state of the interaction field from the ground 
state of the free field. The exact eigenstates can be ex- 
pressed with the "input field" |n_,n+;0) as 



U{0,-T)\n_,n+;0) 
g-i£±oT (^n-,n+ |n_,n+;0)' 



(26) 



6 



for T — > oo (1 + ie). Here, 

U{t,ta)^Texpl^~i dt'H'jit')^ (27) 

is the time-evolution operator in the interaction pic- 
ture, where T is the time-ordering symbol and H'j (t') = 
exp{iHssHt')Hi exj){—iHssHt')- Similarly, we have the 
"output field" 



(n_,n+;0|C/(T,0) 

^-iE±oT {n^, n+; |ri_, n+) 



(28) 



With these formal expressions for the exact eigenstates 
of the total system, we can calculate reduced density ma- 
trix with respect to the AB-subsystem. Then the expec- 
tation value of the projection operators of n+)^^ in 
|n_,7i+) can be calculated as 



Pn 

where 

Q 

and 



_ (n_,n+;0| Q |n_,ri+;0) 
""■^ ~ e-*^±o(2T) |(n_,n+;0 |n_,n+)| 

/ dt'H\ {t' 

J-T 



T i Qn-n+ exp 



(29) 



(30) 



(31) 



Notice that the above equation can not be calculated 
directly, because there is not an explicit expression for 
|n_,n_|-) with respect to the well known basis vectors 
|n_,7i+;0). To get rid of this difficulty, we use the nor- 
malization condition (n_,n-|_ |n_,n^_) = 1 or 



1 = 



(n_,n+;0|£/(T,-T)|n„,n+;0) 
e-iB±o(2T) |(n_,n+;0 |n_,n+)|^' 

Then the Eq. (|29|1 is rewritten as 

p ^ (n^,7i+;0|g|n_,7i+;0) 

(n_,n+;0|[/(r,-r)|n_,n+;0>' 



(32) 



(33) 



If the connection between sites AB and SSH chain is 
switched off, i.e., the hopping integral in Hj^b vanishes, 
states |1, 0; 0) and |0, 1; 0) are the eigenstates correspond- 
ing to vanishing eigen energies. Thus we can expand 
^n_n+, which equals to unity when Hab is exact, as the 
series of g. The Dyson expansion oi U {T, —T) also gives 
the numerator and the denominator in the form of the se- 
ries of g. Comparing the terms of each order with respect 
to small coupling g, the zero order of P„_„^ {Pn}n^) is 

unity, while the first order of Pn_n+ iPn'^n+) is zero as 
a result of the particle number conservation, and second 

order of P, ' ^^^^ 
as 



(Pn_Ti+) is non-zero and can be expressed 



p(2) 

n_n-|- 



^4N 



1 + (n_ — n+) cos k 



I - 1 



-Ok 



(34) 



It approximately describes the difference between the 
eigenstates of Hab and those reduced states from the 
exact eigenstates of the total system. In order to clearly 

(2) 

characterize our model by the correction P„_„^ , we have 

f 2) _____ 

the upper bound of Pn_n+ directly from 1341) as 



p(2) 



< 



(35) 



In the thermodynamic limit, this upper bound becomes 

-,2 



d(2) 



< 



9 



(36) 



This inequality provides us a necessary condition for the 
validity of the effective Hamiltonian Hab ■ Actually, the 
necessary condition is g'^/{2gQ6) <C 1, which requires the 
validity condition 



g/{V2g„) « Vs. 



(37) 



On the other hand, the sufficient condition is g/go <C 1 
for this scheme. Then, if we take S = g/go ^ 1, it 
should satisfy the sufficient and necessary conditions un- 
der which the effective Hamiltonian can work well. 

In this paper, we take S — g/go = 0.01 for the nu- 
merical calculation to demonstrate the validity of this 
scheme. Furthermore, the behavior of Pn_n+ around the 
point (5 = 0, in which region the above condition l|37|) 
is violated, is also investigated by numerical method. It 
shows that Pn_„^ deviate from unity within this region. 
We will discuss this problem in detail in the following. 

Moreover, the mode concurrence between A and B for 
the state |n_,n-|_) can be gained from Pn_n+- It could 
also be used to characterize the validity of the effective 
Hamiltonian Hab- The mode concurrence for the state 
|n_, n+) is 



C„_„^ =max{0,(X)„_„^+2|(r)„_„^|-l}, (38) 

where X = ua + nB — 2nAnB, Y — a^os -I- a^a^i, and 
i^)n^n+ ~ (j^-) 'T'-i- 1 O |rt_, n-|_) denotes the average value 
for an arbitrary operator O Since states \n-,n+) 

and |n_,n+)^^ (g) |0), possess the same symmetry with 
respect to the points A and B, we have 

(^>„_„+ + 2| (i^)„_„, I = (X - 2 (n_ - n+) y)„_„ 
Similarly, we have 



(39) 



Qr 



[X -2 (n_ 



(40) 



Under the validity condition H37|) . P„_„^ is near unity 
and 2Pn_n+ — 1 should be larger than zero. So we have 
mode entanglement 

(41) 



I 


1 - Pio 
(xlO^^) 


p{2) 
10 

(xl0"3) 


1-Poi 
(xlO^^) 


p(2) 
01 

(xlO" 


9 


4.48 


4.53 


5.39 


5.38 


19 


0.41 


0.41 


0.46 


0.47 


29 


4.66 


4.70 


3.69 


3.68 


39 


0.38 


0.38 


0.46 


0.47 


49 


4.60 


4.64 


4.30 


4.30 


59 


0.50 


0.50 


0.45 


0.46 


69 


4.45 


4.49 


5.80 


5.80 


79 


0.67 


0.67 


0.44 


0.44 



TABLE I: Pio and Poi for the finite system with iV = 500, 
Z = lOn + 9, n € [0,7], 5 = 0.01, g = 0.01, and go = 1.0 
obtained by numerical exact diagonalization and from the 
approximately analytical expression Eq. 13411 . They are in 
agreement with each other very well, which shows that the ef- 
fective Hamiltonian Hab can well describe the states of points 
A and B in the ground and first excited states of the total 
system. 



Obviously, when the effective Hamiltonian Hab is valid, 
the sites A and B are maximally entangled for the ground 
and first excited states of the whole system. 

From the Eq. (|21|) about effective hopping integral 
geff, one can see that when I = An + 1 (n = 0, 1, 2, ...), 
|1, 0) is the ground state of the total system, while |0, 1) 
is the ground state for I — An + 3, due to the sign of Qef / 
being positive or negative. We denote the expectation 
value Pn^n+ of the projection operator for the ground 
and first excited states as Pg and Pe, respectively. 

To demonstrate the validity of the effective Hamilto- 
nian of the sites AB, numerical simulation is performed 
for the finite system. We calculate Pg and Pe for the finite 
system with N = 500, I = lOn + 9, n & [0, 7], 6 = 0.01, 
g — 0.01, and = 1.0 by numerical exact diagonaliza- 
tion and from the approximately analytical expression 
Eq. H34|l . The results are listed in Table 1. It shows that 
they are in agreement well and the single-particle (half- 
filled) eigenstates of the effective Hamiltonian Hab can 
well describe the states of points A and B in the ground 
and first excited state of the total system. So Hab is 
valid when temperature T <C AgoS/ks- 

In the above studies, it is found that the magnitude 
of distortion S is crucial for the quantum data bus to 
perform QST and the speed of QST is sensitive to the 
dimerization. There are two types of dimerization cor- 
responding to ± \S\. It is interesting to investigate what 
happens near the transition point (5 = 0. In this region 
the analytical results obtained above is no longer valid 
due to the vanishment of the gap. Numerical simula- 
tion should be employed to calculate geff, Pg and Pe, 
which characterize the property of QST. The numerical 
results obtained by exact diagonalization for the systems 
with N = 100, 200, I = 23, 33, 43 is plotted in Fig. 5. 
It shows that when the systems approach to the critical 
point 6 = 0, the gap AE between the ground and first 
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FIG. 5: The energy gap AE between the ground and first 
excited states, Pg, and Pe of the system with (a) = 100; 
(b) N = 200, and I = 23, 33, 43, obtained by numerical exact 
diagonahzation as the function of S. When S approaches to 
zero, the gap AE has a sharp peak, while Pg tends to 0.25, Pe 
remains unitary approximately. It shows that although QST 
becomes fast, the fidelity decreases rapidly around 5 = 0. 
This observation indicates that the features of two sites A 
and B reveals the transition of the data bus. 



excited states has a sharp peak, while Pg tends to 0.25, 
Pg remains unitary approximately with only a very slight 
jump around the transition point. It indicates that al- 
though QST becomes fast around the critical point, the 
fidelity decreases rapidly. This phenomenon is very simi- 
lar to that discovered by Ref . [23 ■ This fact also implies 
that the features of two probing sites A and B, such as 
the fidelity (mode concurrence) and the recurrent time 
tt/AE, reveal the intrinsic nature of the quantum data 
bus, a quantum phase transition-like behavior of the SSH 
chain. 



QUANTUM DECOHERENCE BY PHONON 
EXCITATIONS 



Finally we consider the decoherence problem due to 
the couplings with auxiliary modes concerning phonon 
excitations. To this end, we need to calculate the ef- 
fect of electron- vibration coupling in the conducting poly- 
mers, which is off-resonate to the electron motion. The 
electron-vibration coupling gives rise to the disadvantage 
of the QST scheme based on SSH chain, though it does 
not dissipate the energy of the electronic subsystem. 

Now we revisit the role of the electron- vibration cou- 
pling Hamiltonian in the process of quantum infor- 
mation transfer based on the virtual excitation of the 
SSH model. We notice that there is no direct interaction 
between the electron motion and the quasi-excitation of 
SSH model when the energy gap is much larger than 
the connection couphngs since [Hab,Ho] = 0. Thus 
the main source of decoherence is due to electron-phonon 
couplings. This will realize a typical quantum decoher- 
ence model for a two level system coupled to a bath of 
harmonic oscillators 0, . 

In order to analyze this decoherence problem more 
quantitatively, we describe the single mode phonon by 
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the perturbation 



(42) 



induced by the quantized vibration of the SSH chain, 
where 7 is a constant and (b) is the bosonic cre- 
ation (annihilation) operator of the phonon excitation. 
In small distortion case, according to H21|l the effective 
Hamiltonian depends on the phonon excitation through 



where 



geff^u;s[l-hib^+b)] (43) 

= (-1)"^ ^. (44) 

Then we can rewrite the total Hamiltonian = Hg + 
Hp + Hq-p describing the quantum decoherence problem 
with 



Hp 

Hq-p 



= UJnb^b, 



(45) 



4uJsJ +b) (Tx:, 



where the Hamiltonian Hq + Hq^p corresponds to Hab 
and has been rewritten in terms of the quasi-spin opera- 
tors 



(46) 



a\^A — a}paB- 



Here, we have ignored the indices of the spin degree of 
freedom for avoiding confusion, because all results we 
obtain in the following are unrelated to spin. 

Obviously, the effective electron-phonon coupling 
Hq-p can only lead to a phase shift in the spin qubit 
initially prepared in the quasi-spin states 



|o,i; 



AB ' 



|i,o) 



AB ■ 



(47) 



It is also pointed out that the ground state (an eigen- 
state of the quasi-spin a^) of the effective Hamiltonian 
is a mode entanglement state. We consider the vibration 
mode initially prepared in a thermal equilibrium state 



|n)(n|. 



(48) 



at the temperature T, where the partition function is 

1 



Z = 



1 



(49) 



and ks is the Boltzman constant. 

Let the electron be initially in a pure state |0) = 
u\^) + v\^). After a straightforward calculation we 
get the density matrix at time t 



p{t) = u{t)m (0i®pp)c/-i(t) 



(50) 



and the corresponding reduced density matrix (t) = 
Trp[p{t)] of AB can be obtained by tracing over the 
phonon modes. The off-diagonal elements of ps (t) can 
be given explicitly as 



p: Wio=Ps it),,=uv*DiT,t)., 



where 



D (T, t) = exp 



2 • 2 / ^ot 
-a sm I — 



(51) 



(52) 



is the so-called decoherence factor P.l2^. The time inde- 
pendent factor 



\ 90(^0 J 



coth 



^0 
2kBT 



(53) 



depends on the temperature T and the distance / between 
A and B. The time dependent behavior of the decoher- 
ence factor is mainly determined by the oscillating factors 
sin^ {ijjQt/2). It indicates that, by increasing I and T, 
becomes larger. It leads to D (T, i) — > except for some 
special instants at 



t 



2m7r 



-,m = 0,1,2. 



(54) 



Here, D {T, t) decays to a minimum Dmin = cxp (— a^) 
as temperature increases. This result shows that thermal 
excitation of the phonons will block a perfect QST to 
some extent and this effect becomes prominent as the 
distance between A and B increases. To consider the 
extent of this influence about the QST scheme based on 
the SSH chain, we can calculate the fidelity or the purity 
P = Trp[ps (t)po (t)] : 



P= l + 2\uv*f [D{T,t) - 1] 



with the ideal state 



-iHABt 



JHABt 



(55) 



(56) 



The Eq. (|55|l shows the intrinsic relation between the 
purity and the decoherence factor. 

Actually the off-diagonal elements are responsible for 
the intrinsic physics of QST, because the time evo- 
lution driven by Hab for A and B will stop if off- 
diagonal elements vanish. We notice that the norms 
of off-diagonal elements is proportional to the decoher- 
ence factor D (T, t) . In order to show the influence of 
the electron- vibration coupling on the QST, 3D diagram 
of D{T,t) is plotted in Fig. 6(a) for the case with 
go = lO^coo, g = lO^coo, 7 = 10"", I = 9, t £ [0,7^^0], 
and T e [fcB/4tJo, lOfcs/t^o]. In Fig. 6(b) and (c), the 
cross sections of the 3D diagram for D (T, t) are plotted 
to demonstrate the behavior of the decoherence factor. 
Actually, there exist other modes of phonon excitations 
in the SSH chain. But in the case of dimerization, the 
mode of phonon we described above dominates the quan- 
tum decoherence and thus it is reasonable that we only 
focus on the dimerization mode. 
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(b) t (C) T 

FIG. 6: (a) A 3D plot of the decoherence factor D (T, t) as 
the function of the temperature and time for the case with 
go = lO^tJo, g = lO^tJo, 7 = 10-^ / = 9, t £ [O.Ttt/wo], and 
T G [ks/^uJo, lOfcs/oJo]- (b), (c) The cross section of D {T,t) 
for different T and I. It shows that thermal excitation of the 
phonons will block a perfect QST as T or Z increases. 



VI. CONCLUSIONS 

In this paper we proposed a QST scheme based on a de- 
sirable quantum data bus, which can be implemented by 



a more practical system, the conducting polymers (poly- 
acetylene) modeled as a SSH chain. This system serves 
as a quantum data bus with the following advantages. 
Firstly, the strength of the induced effective hopping of 
a qubit between two distant sites does not decay rapidly 
over the distance between them, and thus it can realize 
a fast quantum entanglement to transfer quantum infor- 
mation with the always-on coupling and minimal control. 
Secondly, there exists an energy gap to ensure the virtual 
excitation of the electrons in the SSH chain which induces 
the effective hopping of a qubit between two distant sites, 
but does not dissipate their energy. This shows the sim- 
ilar results to that of the spin ladder |0 . Thus a higher 
fidelity can be achieved for our QST scheme based on the 
virtual excitation of the SSH chain in lower temperature. 
Technologically, the more natural and less artificial de- 
signs with practical modulation about coupling constants 
for our protocol push the studies of QST based on real 
physical solid state system. 



Furthermore, based on the analytical and numerical 
results, it is found that the validity of Hab strongly de- 
pends on the distortion 6 of the SSH chain. Thus, as a 
measurable issue, it essentially reflects the intrinsic prop- 
erty of the medium itself. In this sense, such a scheme 
can also be employed to explore the intrinsic property 
of the quantum system. These observations have some 
universality and may motivate us to investigate the func- 
tions of quantum data bus based on other real physical 
systems. 
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